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Abstract 

We present a unique theoretical description of the physics of the spherically trapped iV-atom 
degenerate Fermi gas (DFG) at zero temperature based on an ordinary Schrodinger equation with a 
microscopic, two body interaction potential. With a careful choice of coordinates and a variational 
wavefunction, the many body Schrodinger equation can be accurately described by a linear, one 
dimensional effective Schrodinger equation in a single collective coordinate, the rms radius of the 
gas. Comparisons of the energy, rms radius and peak density of ground state energy are made to 
those predicted by Hartree-Fock (HF). Also the lowest radial excitation frequency (the breathing 
mode frequency) agrees with a sum rule calculation, but deviates from a HF prediction. 
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I. INTRODUCTION 



The realization of the degenerate fermi gas (DFG) in a dilute gas of fermionic atoms 
has triggered widespread interest in the nature of these systems. This achievement com- 
bined with the use of a Feshbach resonance allows for a quantum laboratory in which many 
quantum phenomena can be explored over a wide range of interaction strengths. This leads 
to a large array of complex behaviors including the discovery of highly correlated BCS-like 
pairing for effectively attractive interactions Q^O, ^,0]. While these pairing phenomena 
are extremely interesting, we believe that the physics of a pure degenerate Fermi gas, in 
which the majority of atoms successively fill single particle states, is worthy of further study 
even in the absence of pairing. This topic is addressed in the present paper, which follows 
a less complete archived study of the topic. [f| 

The starting point for this study is that of a hyperspherical treatment of the problem 
in which the gas is described by a set of 3N — 1 angular coordinates on the surface of a 
3iV dimensional hypersphere of radius R. Here N is the number of atoms in the system. 
This formulation is inspired by a similar study of the Bose-Einstein condensate, p, 0, 0| 
Furthermore these same coordinates have been applied to finite nuclei. [n| The formulation 
is a rigorous variational treatment of a many-body Hamiltonian, aside from the limitations 
of the assumption of pairwise, zero-range interactions. In this paper, we consider only 
filled energy shells of atoms. This is done for analytic and calculational simplicity, but the 
treatment should apply to any number of atoms in open shells with modest extensions that 
are discussed briefly. 

The main goal of this study is to describe the motion of the gas in a single collective 
coordinate R, which describes the overall extent of the gas. The benefit of this strategy 
is that the behavior of the gas is reduced to a single one-dimensional linear Schrodinger 
equation with an effective hyperradial potential. The use of a real potential then lends itself 
to the intuitive understanding of normal Schrodinger quantum mechanics. This method also 
allows for the calculation of physical quantities such as the energy and rms radius of the 
ground state; these observables agree quantitatively with those computed using Hartree- 
Fock methods. The method also yields a visceral understanding of a low energy collective 
oscillation of the g breathing mode. 

Beyond the intuitive benefits of reducing the problem to an effective one-dimensional 
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Schrodinger equation, another motivation for developing this hyperspherical viewpoint is 
that it has proven effective in other contexts for describing processes involving fragmen- 



tation or collisions in few-body and many-body systems. [111 . Il2j Such processes would be 
challenging to formulate using field theory or RPA or configuration interaction viewpoints, 
but they emerge naturally and intuitively, once the techniques for computing the hyper- 
spherical potential curves for such systems are adequately developed. To this end we view 
it as a first essential step, in the development of a more comprehensive theory, to calculate 
the ground-state and low-lying excited state properties within this framework. Then we 
can ascertain whether the hyperspherical formulation is capable of reproducing the key re- 
sults of other, more conventional descriptions, which start instead from a mean-field theory 
perspective. 

The paper is organized as follows: Section II develops the formulation that yields a 
hyperradial ID effective Hamiltonian; in Section III we apply this formalism to a zero-range 
s-wave interaction and find the effective Hamiltonian in both the finite N and large N 
limits. In Sections Ilia and Illb we examine the nature of the resulting effective potential 
for a wide range of interaction strengths and give comparisons with other known methods, 
mainly Hartree-Fock (HF); Section IIIc is a brief discussion of the simplifications that can 
be made in the limit where N — > oo; finally in Section IV we summarize the results and 
discuss future avenues of study. 



II. FORMULATION 

The formalism is similar to that of reference [jj, but we will reiterate it for clarity and 
to make this article self-contained. Consider a collection of N identical fermionic atoms 
of mass m in a spherically symmetric trap with oscillator frequency u), distributed equally 
between two internal spin substates. The governing Hamiltonian is 

# = ^E V *+ E * + E u ** ft) (2- 1 ) 

i=l i=l i>j 

where Ui nt if) is an arbitrary two-body interaction potential and — fi — r}. We ignore 
interaction terms involving three or more bodies. In general, the Schrodinger equation 
that comes from this Hamiltonian is very difficult to solve. Our goal is to simplify the 
system by describing its behavior in terms of a single collective coordinate. To achieve this 
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aim we transform this Hamiltonian into a set of collective hyperspherical coordinates; the 
hyperradius, R, of this set is given by the root mean square distance of the atoms from the 
center of the trap. 

Rs (ji^) ■ < 22) 

So far we only have one coordinate for the system, but we need to account for all 3N spatial 
coordinates and also the spin degrees of freedom. This leaves 3iV — 1 angular coordinates 
left to define. We have 2N of the angles as the independent particle spherical polar coordi- 
nate angles (9i , <j><\ , 9p., 0a, ...,9 N ,(p N ). The remaining N — 1 hyperangles are chosen by the 
convention used in [10], which describes correlated motions in the radial distances of the 
atoms from the trap center, 



V 1 r 2 
Z-/j=i j 

tanctj = — , (2.3) 

r i+ i 

i = 1,2,3, ...,N- 1. 

Alternatively we may write this as 

N-l 

r n = \f~NR cos a n -i Y\ sina^ (2.4) 

j=n 

0<«i<|, j = l,2,...,iV-l 

N-l 

where we define cos«o = 1 an d ]~[ sin<x,- = 1. For the purposes of this study, the set 

j=N 

of 3iV — 1 hyperangles will be referred to collectively as Q. The particular definition of 
the hyperangles will not play a significant role in the actual formalism, but we give the 
definitions here for completeness. 

After carryin g ou t this coordinate transformation on the sum of Laplacians, the kinetic 
energy becomes [13j 



2M 



I f R 3N-l 9 \ A 



,2 



(2.5) 



R^^dRy dRJ R 2 
Here M = Nm and HA is the grand angular momentum operator which is similar to 
the conventional angular momentum operator and is defined by 

a2 = - E 4> ^ ■ r, Ji7~ x >Ji7 : (2 - 6) 

i>j 3 1 
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for all cartesian components, Xi, of the 3N dimensional space. The isotropic spherical 
oscillator potential becomes simply 



i A i 



-^Y- = -Mu 2 R 2 . (2.7) 



2 

i=l 



The sum of Eqs. 12.51 and 12.71 gives a time- independent Schrodinger equation, (R, ft) 
E<Z> (R, ft), of the form 







-V ( d 2 , (3N - 1) (3N - 3) A 2 ^ , 1 



2M V &R 2 2# 2 R 2 2 



+ " ^ " - « + (2-8) 



+ ^ ^int (r* - r}) - 



where \1/ (i?, ft) has been multiplied by i?( 3Ar_1 )/ 2 to remove first derivative terms in the 
hyperradius. These mathematical transformations have not yet accomplished much. We 
started with a 3N dimensional Schrodinger equation to solve, and we still have a 3N di- 
mensional Schrodinger equation. To simplify, we assume that \I/ is approximately separable 
into an eigenfunction of the operator A 2 , a "hyperspherical harmonic" . multiplied by an 
unknown hyperradial function. Hyperspherical harmonics (HHs, see LJ for more details) 
are generally expressed as products of Jacobi polynomials for any number of dimensions. 
Their eigenvalue equation is 

A 2 $ Ai , (ft) = A (+A + 3N - 1) $ A „ (ft) (2.9) 

Where A = 0, 1, 2, ... and v (omitted below for brevity) stands for the 3N — 2 other quantum 
numbers that are needed to distinguish between the (usually quite large) degeneracies for a 
given A. The separability ansatz implies that 

* (R, ft) = F (R) $ A (ft) . (2.10) 

Here we assume that $ A (ft) is the HH that corresponds to the lowest value of the hyper- 
angular momentum that is allowed for the given symmetry of the problem, i.e. that is 
antisymmetric with respect to interchange of indistinguishable fermions. This choice of trial 
wavefunction indicates that we expect the overall energetics of the gas to be described by its 
size; as such, fixing the hyperangular behavior is equivalent with fixing the configuration of 
the atoms in the gas. In nuclear physics this is known as the K harmonic method. Alterna- 
tively, this may be viewed as choosing a trial wavefunction whose hyperradial behavior will 



be variationally optimized, but whose hyperangular behavior is that of a non-interacting 
trap dominated gas of fermions. 

To utilize this as a trial wavefunction, we evaluate the expectation value ($a \H\ $a), 
where the integration is taken over all hyperangles at a fixed hyperradius. This approach 
gives a new effective linear ID Schrddinger equation H e ffR^ 3N ~ 1 ^ 2 F (R) = ER^'^^F (R) 
in terms of an effective Hamiltonian H e ff given by 

-h 2 ( d 2 K(K+1)\ 1 2 2 . . v 

Here AT = A + 3 (N - 1) /2. 

We now must force our trial wavefunction to obey the antisymmetry condition of fermionic 
atoms. To antisymmetrize the total wave function F (R) $a note that Eq. 12.21 indicates 
that R is completely symmetric under all particle coordinate exchange, thus the antisym- 
metrization of the wavefunction must only affect <3>a • Finding a completely antisymmetric 
<3>a (^) for any given A is generally quite difficult and is often done using recursive techniques 



s) or using 
191 We use 



like coefficient of fractional parentage expansions (see lfl [lj, [l5( for more detai 
a basis of Slater determinants of independent particle wave functions. 12, 3, 
a simplified version of the second method combined with the following theorem, proved in 
reference [3] and developed in Appendix B. 

Theorem 1 The ground state of any non-interacting set of N particles in an isotropic 
oscillator is an eigenfunction of A 2 with minimal eigenvalue A (A + 3N — 2) where A is 
given by the total number of oscillator quanta in the non-interacting system. 

where E NI is the total ground state energy of the non-interacting N-body system. 

With Theorem 1 in hand we may find $a (^) in terms of the independent particle coor- 
dinates. The non-interacting ground state is given by a Slater determinant of single particle 
solutions. 

N 

Vni (n,r 2 , ...,r N ,a 1 ,a 2 , -,crjv) = (- 1 ) PjP II jRn ^ ( r i)yi imi \rn s .) . (2.13) 

P 8=1 

Here o~i is the spin coordinate for the ith atom, R ni ii is the radial solution to 
the independent particle harmonic oscillator for the ith particle given by rR n £ (r) = 



N n £exp (— r 2 /2/ 2 ) (r/l) i+1 Li +1//2 [(r/7) 2 ] where L® (r) is an associated Laguerre polynomial 
with I = ^Jh/muj. ya imi (u) is an ordinary 3D spherical harmonic with uoi as the spatial solid 
angle for the zth particle, \m s .) is a spin ket that will allow for two spin species of atoms, 
||) and ||). The sum in Eq. 12.131 runs over all possible permutations P of the iV spatial 
and spin coordinates in the product wavefunction. We now apply Theorem 1 which directly 
leads to a ty^i that is separable into a hyperangular piece and a hyperradial piece 

^ni (ri, r 2 , f N , oi, a 2 , ...a N ) = G (R) $ A (Q, a u a 2 , ...cr N ) . (2.14) 

Here G (R) is a the nodeless hyperradial wave function describing the ground state of N non- 
interacting atoms in an isotropic oscillator trap. A derivation of G (R) is given in Appendix 
A: 

R (3N-1)/2 G{R) = Axexp (_ R 2/ 2C *} ( | j , (2.15) 

here A\ is a normalization constant and C = yJh/Mu = l/y/N. Combining Eq. 12.141 with 
Eq. 12. 151 now gives us $a (^) 

N 

E {~ l T P 11 Rrnti in) yttm I™*) 
$ A (fi, oi, 2 , ...0-jv) = i=1 (2.16) 

where we must make the variable substitutions in the numerator using Eqs. 12.21 and 12.31 In 
the following, for brevity, the spin coordinates (o"i, ...,<tn) will be suppressed, i.e. $ A (0) = 
$a (^, 02 > ■••0"jv). Interestingly, this must be a function only of the hyperangles and thus 
all of the hyperradial dependence must cancel out in the right hand side of 12.161 We are 
now ready to start calculating the interaction matrix element ($a Wi n t \ 3>a}- 

III. ZERO-RANGE S-WAVE INTERACTION. 

Here we specify Ui n % (r) as a zero-range two body interaction with an interaction strength 
given by a constant parameter g. 

U mt (r)=g5 3 (r) (3.1) 

For a small two body, s-wave scattering length a we know that g = .I20I For stronger 

m 

interactions, i.e. \kta\ > 1, this approximation no longer holds and g must be renormalized. 
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Now we must calculate the effective interaction matrix element given by 

U eff (R) = gJ2(®x\t 3 (r il )\$x)- (3.2) 

i>j 

Degenerate ground states cause some complications for this formulation which we avoid 
by restricting ourselves to the non-degenerate ground states that correspond to filled energy 
shells of the oscillator, i.e. "magic numbers" of particles. With moderate extensions the 
degeneracies can be taken into account by creating an interaction matrix, but the magic 
number restriction should still give a good description of the general behavior of systems 
with large numbers of atoms. The total number of atoms and the hyperangular momentum 
quantum number A are most conveniently expressed in terms of the number n of single 
particle orbital energies filled: 

N = " ( " +1 > ( " + 2) (3.3a) 
X = '"- 1 >"'" 4 +1 "" + 2 » (3.3b) 

where kf is the peak non-interacting Fermi wave number. In the limit where iV 3> 1, we 
write A and kf in terms of the total number of particles N, 

x ( 3iV ) 4/3 

A -> (3.4a) 

*/ - ^ (3iV) 1/3 (3.4b) 

Next we combine I3.4al with $^ (^) from !2.1o1 and calculate the interaction matrix element. 

Since $a (^) is antisymmetric under particle exchange we may do a coordinate transpo- 
sition in the sum r*j — > f 2 and fj — > fx. Each transposition pulls out a negative sign from <3>^ 
and we are left with 



U eff (R)=gJ2(®x\S 3 (^2i)|$A> 



i>j 



. g ^i) (4A |, 3(?2l) | $A> . 

Appendix C details the calculation of the matrix element ($a (^*2i)| ^a)- The result is 

U *ff ( R ) = aj^ks (3-5) 
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where Cat is a constant that is dependent only on the number of atoms in the system. 
While Cn has some complex behavior for smaller numbers of particles, we have seen that 
for iV > 100, CV quickly converges to 

^2 32V 7 / 2 / 0.049 0.277 



C 



N 



3 35tt 3 
0.02408iV 7/2 



1 + 



AT2/3 

0.049 



AT4/3 

0.277 



+ ... 



+ ... 



(3.6) 



jV 2 / 3 V 4 / 3 

where the higher order terms in 1/N were found by fitting a curve to numerically calculated 
data. Values of C N /N 7 ^ ' - - Fig. □ shows both the 




0.0005 0.001 0.0015 0.002 0.0025 
1/N 



FIG. 1: Values of C N divided by the large N limit, C N 0.02408V 7 / 2 , versus 1/N are shown. 
The circles are the calculated value while the curve is the fit stated in Eq. 



calculated values of and the values from the fit in Eq. 13.61 versus 1/N. In both the table 
and the plot we can clearly see the convergence to the large V value of C N « 0.02408V 7/2 . 

We now may write the effective one dimensional, hyperradial Schrodinger equation, 
if e// i? (3Ar - 1)/2 F (R) = ERW-WF (R) where H eff is given by 

*«" = -S?5F + ir «'W < 37 > 

where V e ff (R) is an effective hyperradial potential given by 

We reinforce the idea that this effective Hamiltonian describes the fully correlated motion of 
all of the atoms in the trap albeit within the aforementioned approximations. As expected, 
for V = 1, with C\ = 0, Eq. 13.81 reduces to a single particle in a trap with K as the angular 



n 


N 


A 


C N /N 7 / 2 


1 


2 





5-5 ~ 0.0127 


2 


8 


6 


0.0637 


3 


20 


30 


0.0251 


4 


40 


90 


0.0244 


5 


70 


210 


0.02435 


15 


1360 


14280 


0.0241 


30 


9920 


215760 


0.02409 


100 


343400 


25497450 


0.02408 



TABLE I: N, A and C N /N 7 / 2 for the several filled shells. We can see that C N /N 7 / 2 quickly 
converges to the Thomas-Fermi limiting value of 32y / 2/3/357r 3 ~ 0.02408 to several digits. 

momentum quantum number £. Note that the form of V e ff is very similar to the effective 
potential found for bosons by the authors of |7(. What may be surprising is the extra 
term of 3 (N — 1) /2 contained in K. The kinetic energy term in V e ff is controlled by the 
hyperangular momentum, which in turn reflects the total nodal structure of the iV-fermion 
wavefunction. This added piece of hyperangular momentum summarizes the energy cost of 
confining N fermions in the trap. This repulsive barrier stabilizes the gas against collapse 
for attractive interactions, i.e. g < 0. 

A final transformation simplifies the radial Schrodinger equation, namely setting E = 
E NI E' and R = ^J(R 2 ) NI R' where 

E NI = hu + -rpj > (3.9a) 

are the non-interacting expectation values of the energy and squared hyperradius in the 
ground state. In the following any hyperradius with a prime, R' denotes the hyperradius in 
units of a/ (R 2 ) ni . Under this transformation, with the use of Eqs. I3.4al and I3.4b[ and in 
the limit where — > oo, the hyperradial Schrodinger equation becomes 
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where m* = (A + 3N/2) 2 . The effective potential takes on the simple form 



Veff (R') 



E 



1 . l R a , ° k n 

2R' 2 2 i?' 3 ' 



(3.11) 



where a = 1024/28357T 3 and 7 = g/hul 2 . We note that the only parameter that remains in 
this potential is the dimensionless quantity fc/7 and that in oscillator units 7 = g. 
In the limit where iV — > 00 the effective mass m* becomes 



m 



J_ {3N f* 



(3.12) 



For large numbers of particles, the second derivative terms in Eq. 13. 101 becomes negligible, 
a fact used later in Section IIIc. 

The behavior of V e ff (R') versus R! is illustrated in Fig. El for various values of fc/7. For 




0.5 1 1.5 2 2.5 

R/V < R 2 >ni 



FIG. 2: The effective potential as a function of the hyperradius for kfj = (solid), kfj = 15 
(dashed) and kfj = —5. (dot-dash) 

kf'y = (solid curve), the non-interacting limit, the curve is exact and the ground state 
solution is given by equation 12.151 For non-zero values of g, V e ff acquires an attractive 
(/c/7 < 0) or repulsive (kfy > 0) 1/R' 3 contribution as indicated by the dot-dashed and 
dashed lines respectively. For k^ < the DFG is metastable in a region which has a 
repulsive barrier which it may tunnel through and emerges in the region of small R' where 
the interaction term is dominant. It should be noted, though, that small R' means the overall 
size of the gas is small. Thus the region of collapse corresponds to a very high density in 
the gas. In this region several of the assumptions made can fall apart, most notably the 



assumption dealing with the validity of the two- body, zero-range potential. 



For kf'y > 



the positive 1/i?' 3 serves to strengthens the repulsive barrier and pushes the gas further out. 
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A. Repulsive interactions (g > 0) 



For positive values of the interaction parameter g we expect the predicted energy for this 
K harmonic method to deviate from experimental values, since the trial wavefunction does 
not allow any fermions to combine into molecular pairs as has been seen in experiments. |l], 0, 
0, |^ Our method only can describe the normal degenerate fermi gas. The strong repulsive 
barrier for repulsive interactions shown in Fig. |21 arises as the gas pushes against itself 
which increases th( Q] and 0] compare 



1.15 



E9 



1.05 




FIG. 3: The ground state energy in units of the non-interacting energy versus kfy for 240 atoms 
calculated using the K harmonic method (curve) and using Hartree-Fock (circles). Inset: the 
difference in the ground state energies predicted be the K harmonic (Ek) and Hartree-Fock (Ehf)- 
Clearly the K harmonic energies are slightly higher than Hartree-Fock. 

the ground state energy and average radius squared respectively of 240 trapped atoms, 
plotted as a function of kfj with a Hartree-Fock (HF) calculation. The inset in Fig. El 
shows that the K harmonic energies are slightly above the HF energies; since both methods 
are variational upper bounds, we can conclude that the Hartree-Fock solution is a slightly 
better representation of the true solution to the full Schrodinger equation with ^-function 
interactions. 

An added benefit of the K harmonic method is that we now have an intuitively simple 
way to understand the energy of the lowest radial excitation of the gas, i.e. the breathing 
mode frequency. Fig. El shows that as kfy increases the repulsion increases the curvature 
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5 10 15 



FIG. 4: The ground state average squared radius of the gas atoms in units of the non- interacting 
rms squared radius is plotted versus kfj. The calculations considered 240 atoms in both the K 
harmonic method (curve) 




I i I i I i I i L 

2 4 6 8 
k f 7 



FIG. 5: The lowest breathing mode excitation (ljq) in units of the trap frequency is plotted versus 
kfj for the K harmonic method (solid curve) and for the sum rule (circles). Also shown as 
dashed curves are the lowest eight radial excitation frequencies predicted in the Hartree-Fock 
approximation . 

at the local minimum, whereby stronger repulsion causes the breathing mode frequency to 
increase. Fig. |5] compares the breathing mode frequency calculated using the K harmonic 
method to the sum rule prediction '22] based on HF orbitals, and also the lowest eight radial 
excitation frequencies predicted by Hartree-Fock. As anticipated, the K harmonic method 
and the sum rule method agree that the breathing mode frequency will increase with added 
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repulsion. Interestingly both the K harmonic and sum rule methods disagree qualitatively 
with all eight of the lowest HF excitations. This difference is attributed to the fact that 
Hartree-Fock on its own can only describe single particle excitations while both the sum rule 
and the K harmonic methods describe collective excitations in which the entire gas oscillates 
coherently. 

As another test of the K harmonic method we calculate the peak density of the gas. To 
do this we first define the density 

Pif) = j (ll^^ 3 (f 4 -r^| 2 . (3.13) 

It can be seen that integration over r using this definition gives j p (r) d 3 r = N. We recall 
that our separable approximation takes the form \1/ = F (R) $a (^)- Use of 13.131 and the 
antisymmetry of $^ gives 

p(Q)=N J dRR 3 "- 1 \F{R)\ 2 J dtt5 3 (r N ) |$ A (fi)| 2 . (3.14) 

The hyperangular integration is carried out in Appendix D. The result is 

" 0) -J^« (3.15) 



Pni (0) J R 3 

/ 3 p f\ _l_ SN/2) 

where p^i (0) is the non-interacting peak density and £ = , — . For 

N d ' 2 T (A + 3 (N — 1) /2) 

the nth filled energy shell the non-interacting peak density is given by 



PNI (0) 



A f 



f 



2mw 



n + l/2) 



3/2 



6tt 2 6ty z 

Note that the peak density is not given by R 3N ~ l \F (R)\ z evaluated at R = 0, this describes 
the probability of all of the particles being at the center at once. 

Fig. El compares the K harmonic and HF peak densities. Clearly the density does decrease 
from the non-interacting value. The two methods are in good qualitative agreement, but 
Hartree-Fock seems to predict a slightly lower density, presumably a manifestation of the 
slightly inferior K harmonic wavef unction, F (R) <& \ (Q) . 



B. Attractive interactions (g < 0) 

In this section we examine the behavior of the gas under the influence of attractive s-wave 
interactions < 0) . For attractive interactions the gas lives in a metastable region and 
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FIG. 6: The peak density in units of {h/mui) ' versus kfy predicted by the K harmonic (solid) 
and HF (circles). Both sets of calculations were done for a filled shell of 240 atoms. 




0.5 1 1.5 2 

R/V < R 2 >ni 



FIG. 7: V e ff versus R' for several values of kfy. kfy = kfy c (dashed) and from top to bottom 
kfy = -5.3, -8.3, -11.3, -19.3 (all solid) 

can tunnel through the barrier shown in Fig. |21 Fig. [7| shows the behavior of V e ff for 
several values of kfy. The location of the local minimum gets pulled down with stronger 
attraction as the gas pulls in on itself and deeper into the center of the trap. Further, as 
the strength of the interaction increases, the height of the barrier decreases. In fact beyond 
a critical interaction strength 7 C the interaction becomes so strong that it always dominates 
over the repulsive kinetic term. At this critical point the local extrema disappear entirely 
and the gas is free to fall into the inner "collapse" region. The value of 7 C can be calculated 
approximately by finding the point where V e ff loses its local minimum and becomes entirely 
attractive. This is not exact as the gas will have some small zero point energy that will 
allow it to tunnel through or spill over the barrier before the minimum entirely disappears. 
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FIG. 8: The ground state energy (in units of the non-interacting energy) versus kfj for 240 atoms 
calculated using the K harmonic (curve) and Hartree-Fock (circles) methods. Inset: The difference 
in the ground state energies predicted be the K harmonic (Ek) and Hartree-Fock (Ehf)- Clearly 
the K harmonic prediction is slightly higher. 

This critical interaction strength is given by 

1897T 3 1 

kfic = 777 ~ —15.31. 

} ' 256 5 x /4 

Just before the minimum disappears, its location is given by -R^in = 5 -1 / 4 , with an energy 
of V e ff (R' min ) = VbE NI /3 ~ 0.75E NI . This means that if the gas is mechanically stable 
for all values of the two body scattering length, i.e. a — > — oo, in this approximation, there 
must be a renormalization cutoff in the strength of the 5-function such that fe/7 > —15.31 
for all a. With this in mind, we begin to examine the behavior of the DFG for the allowed 
values of kfj. 

Figs. |H1 and El show a comparison of the ground state energy and rms radius of the gas 
versus kfj down to the k^ c as calculated in the K harmonic and Hartree-Fock methods. 
Again Hartree-Fock does just slightly better in energy, which we interpret as Hartree-Fock 
giving a slightly better representation of the actual ground state wavefunction. The energy 
difference becomes largest as the interaction strength approaches the critical value. This 
increase is due to the fact that Hartree-Fock predicts that collapse occurs slightly earlier 
with kf'-fc « — 14.1. As the interaction strength increases, the energy and rms radius of 
the gas decrease. As kfy approaches kfj c the overall size of the gas decreases sharply and 
from Eq. 13.151 we expect to see a very sharp rise in the peak density. This behavior is 
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FIG. 9: The average squared radius of the Fermi gas ground state in units of the non- interacting 
value is plotted versus kfy. The calculations are for 240 atoms in both the K harmonic method 
(curve) and Hartree-Fock 




k f a 

FIG. 10: The peak density in units of (h/muj)^ 3 ^ 2 versus kfy predicted by the K harmonic (solid 
line) and HF (circles). Both calculations were carried out for a filled shell of 240 atoms. 

apparent in Fig. which displays the peak density versus kfy for both the K harmonic 
and Hartree-Fock. 

While the local minimum present in V e ff only supports metastable states, it is still 
informative to examine the behavior of the energy spectrum versus fc/7, beginning with 
the breathing mode frequency. As the interaction strength becomes more negative Fig. 
shows that the curvature about the local minimum in V e ff decreases. This "softening" of the 
hyperradial potential leads to a decrease in the breathing mode frequency in the outer well. 
Fig. [TT1 shows the breathing mode vs. kfy predicted by the K harmonic (curve) method 
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FIG. 11: The frequency of the lowest energy radial transition in units of the trap frequency versus 
kfy predicted by the K harmonic method (solid line) and by the sum rule (circles). Also shown 
are the lowest eight radial transitions predicted by Hartree-Fock. 

and also using the sum rule with Hartree-Fock orbitals (circles). Also shown in Fig. HP are 
the lowest eight Hartree-Fock excitation frequencies for a filled shell of 240 atoms. Again, 
the K harmonic method agrees quite well with the sum rule, while both differ qualitatively 
from the HF prediction. The sharp decrease in the breathing mode frequency that occurs as 
kfj — ► fc/7 c is a result of the excited mode "falling" over the barrier into the collapse region 
as the barrier is pulled down by the interaction. Figure H~2l displays some energy levels in the 
metastable region as functions of kfy near kfj c . Because of the singular nature of the 1/R 3 
behavior in the inner region we have added an inner repulsive 1/R 12 barrier to truncate 
the infinitely many nodes of the wavefunction in the inner region. The behavior of the 
wavefunction is not correct within this region anyway because recombination would become 
dominant, and in any case the zero-range interaction is suspect beyond \kfa\ ~ 1 and it 
must be renormalized. Figure shows three distinct types of energy level. Levels that are 
contained in the local minimum (shown in blue) are decreasing, but not as quickly as the 
others; levels that are in the collapse region (shown in red) have a very steep slope as they 
are drawn further into toward R = 0; and energy levels that are above the barrier in V e ff 
(shown in purple) have wavefunctions are in both the collapse region and the local minimum. 
As kfy decreases, the higher energy levels fall over the barrier into the collapse region earlier, 
until finally just before kfy c is reached the lowest metastable level falls below the "ground 
state". This corresponds to the breathing mode behavior seen in Fig. Of course all 
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FIG. 12: (color online) A portion of the energy spectrum vs kfy close to the critical point = 
kfj c . Levels in the metastable region (blue) decrease slowly while levels in the collapse region (red) 
decrease very quickly. Energy levels above the barrier in V*tt (purple) live both in the collapse 
region and the metastable region. 

of this applies only if there is no further hyperradial dependence of g. If the interaction 
becomes density dependent, as is the case in some types of renormalization, a change in 
the hyperradius will change the density and thus the interaction coupling parameter, i.e. 



C. The large N limit 

In Sections Ilia and Illb the ground state energy and expectation values discussed were 
found by solving the hyperradial effective Schrodinger equation 13.71 for a finite number of 
particles. Here we discuss the behavior of the N fermion system in the limit where N is 
large. To do this we exploit the fact that the d 2 /dR 2 term in the effective Hamiltonian, 
the "hyperradial kinetic energy", becomes negligible. In this limit we see that the total 
energy of the system is merely given by E = V e ff (-Rmm), as is the case in dimensional 
perturbation theory. |23| To find the ground state energy we must merely find the minimum 

dV ef f 

(local minimum for g < 0) value of V e ff. Accordingly, we find the roots of j; = 0. Using 

dR' 

Eq. 13.111 we simplify this to 

kfl = ^-R' mi n {Kin ~ 1) (3-16) 

where R' min is the hyperradial value that minimizes V e ff. The solutions to Eq. 13.161 are 
illustrated graphically in Fig. HH1 for any given kfy we need only look for the value of R' min 
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FIG. 13: Plot of the potential curve extrema which obey (|3.16j) . Examination of this plot tells us 
the behavior of R'^^ for all allowed values of kfy including the existence of the critical point kff c 
located at the minimum of the plot where the maximum and minimum coincide. 

that gives that value. The exact solution of Eq. 13.161 cannot be determined analytically 
for all values of k/j, but Fig. fTTTl shows that for kfy > there is always only one positive, 
real R' min that satisfies Eq. 13.161 This corresponds to the global minimum discussed for 
repulsive interactions. For kfy < things are a bit more complicated. Fig. IT31 shows that 
for kfjc < < there are two solutions to 13.161 The inner solution is a local maximum 
and corresponds to the peak of the barrier seen in Fig. [7| the outer solution corresponds 
to the local minimum where the DFG lives. The local minimum is the state that we are 
concerned with here as this will give the energy and hyperradial expectation values of the 
metastable Fermi gas. The value of kfy where these two branches merge is the place where 
the local maximum merges with the local minimum namely the critical value fc/7 c . Any 
value of less than kfy c has no solution to 13.161 and thus we cannot say that there is a 
region of stability. For kfy = 0, the non- interacting limit, we see that there are two solutions 
R' min = and R' min = 1. The solution -R^in = must be discounted as there is a singularity 
in V e ff at R = 0. Thus in the non-interacting limit R' min — > 1 and V e ff (i? m i n ) — > E NI , as 
expected. 

Substitution of 13. 161 into 13 . 1 II gives the energy of the ground state in the large iV limit, 
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FIG. 14: The percentage difference between the energy found by minimizing V e ff and the energy 
found by explicitly solving the hyperradial Schrodinger equation for 240 atoms. 

as a function of the size of the gas. 

V eff (i? min ) _ 1 + 5R\ 



/4 

min 



En j 6R' 2 



111111 



These solutions to Eq. 13.161 immediately give the ground state energy of the gas versus 
kfy. Fig. El shows the percentage difference of the ground state energy found by this 
minimization procedure and that of 240 particles found by diagonalizing H e ff in Eq. 13.71 
We see in Fig. H"4l that for \kfy\ < 10 the energy difference is less than 0.5%. As |fc/7| 
gets larger the difference increases, but in the range shown the difference is always less than 
±1.5%. 

Another result from l3.l01 is the fact that in the large N limit the commutator [H e ff, R] — > 

0. Thus for any operator that is solely a function of the hyperradius O (R) the ground 
state expectation value in the large N limit is given by the operator evaluated at R' min , 

1. e. (O (R')\ = 0(R' min ). This tells us that the large N limit wavefunction is given by 
[i?( 3Ar_1 ^ 2 G (i?)] 2 = 5 (R — Rmin) • We can perturb this slightly and say that the ground 
state hyperradial wavefunction can be approximated by a very narrow Gaussian centered 
at -R m in- To find the width of this gaussian we approximate K// about -R m i n as a harmonic 
oscillator with mass m* and frequency ujq. We may find ujq by comparing the oscillator 
potential with the second order Taylor series about R m i n in V e ff , i.e. 



i 



1 1 / d 2 V f 



eff 



m* E NI \ dR' 2 
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(3.17) 



R'=R' - 

mm / 



The effective hyperradial oscillator length of R'^ N ~ X ^ 2 G (R) is then given by Iq = ^h/m*uo. 

The breathing mode frequency is now simply found to be ujq. The frequency in Eq. 13.171 
is in units of the non-interacting energy, so to get back to conventional units we multiply by 
E NI /h. From Eq. ETTUl we know that m* = mE NI N (R 2 ) m /H 2 . Noting that N(R 2 ) NI = 
l 2 E]s[i/hjJ this leads to 



d 2 V f 



(3.18) 



R'=R' - 

mm y 



E NI \ OR' 2 

for Uq, the breathing mode frequency in units of the trap frequency. Using Eq. 13.111 and 
substituting in Eq. !3.16l to evaluate at the minimum gives that 



R 

min 



We note that this is now dependent only on the value of kff, i.e. for a fixed kf'j the predicted 
breathing mode is independent of the number of atoms in the system in the large N limit. 
This prediction can be compared with that predicted by the sum rule using the formula 



found by the authors of Ref. 
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^ = a/4+^- (3.19) 



Here E int and E^ are the expectation values of the interaction potential and the oscillator 
potential in the ground state respectively. If we insert the expectation values predicted here 
by the K harmonic method we find 



R ° 

min 



Substituting 13.161 for fc/7 gives 



R 



linn 



which agrees exactly with the frequency predicted in Eq. 13.191 Fig. ITB1 shows the breathing 
mode frequency predicted by Eq. 13. 191 versus kf^y. We see the same behavior in this plot as 
was seen in Fig. ITT1 where the breathing mode frequency dives to zero as — > fc/7 c - 



IV. SUMMARY AND PROSPECTS 



We have demonstrated an alternative approach to describing the physics of a trapped 
degenerate Fermi gas from the point of view of an ordinary linear Schrodinger equation 
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FIG. 15: The breathing mode (in oscillator units) in the large N limit versus kfj. Note that as 
kfj — » kfj c the frequency drops to zero as the local minimum disappears. 

with two body, microscopic interactions. The use of a hyper spherical variational trial wave 
function whose hyperangular behavior is frozen to be that of the K harmonic yields a one- 
dimensional effective potential ()3.8)1 in a collective coordinate, the hyperradius R. This 
approach yields an intuitive understanding of the energy and size of the DFG in terms of fa- 
miliar Schrodinger quantum mechanics. Perhaps surprisingly, this approximation gives good 
results for the ground state energy, rms radius and peak number density of the gas. These 
are all in quantitative agreement with HF, while the breathing mode frequency compares 
well with that found using the sum rule method, but gives a qualitative improvement over 
the lowest calculated HF radial excitation frequencies. 

The work presented here has been limited to the case of filled energy shells, i.e. "magic 
numbers" of atoms, and it should be readably generalizable, with minor extensions, to open 
shells as well. The present results are limited to a spherically symmetric trap. Generalization 
to an cylindrically symmetric "cigar" trap should be possible with a judicious choice of 
coordinates, and will be presented elsewhere. 

For strongly attractive interactions, this picture predicts an instability in the DFG in 
a manner qualitatively similar to the physics of the "Bosenova" . j2j For the gas to remain 
stable across the BEC-BCS crossover regime the strength of interaction potentia 
bounded from below, fc/7 > —15.31. Preliminary results from another study 
that renormalization of the singular 5-function interaction accomplishes precisely that, and 
apparently prevents collapse. The full interrelation between this picture and that of pairing 



must be 



24( indicate 
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in the BEC-BCS crossover region is beyond the scope of this study and is a subject that 
will be relegated to future publications. 
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APPENDIX A: HYPERRADIAL SOLUTION TO THE NON-INTERACTING 
OSCILLATOR 

We wish to derive Eq. 12.151 the nodeless hyperradial solution to the Schrodinger equation 
for N non-interacting particles in an isotropic oscillator. The Schrodinger equation for this 
system is given by 

T N /_fi2 1 \ 

Et v « + H - £ 117 = (A1) 

where r 1; r 2 , are the cartesian coordinates for each atom from the trap center. 

To begin we examine the radial Schrodinger equation for a single particle in an isotropic 
trap 

-h 2 ( d 2 £(£ + l)\ 1 22 \ , , . n , AO , 
' 1 + -mu 2 r 2 - E ni rf ne (r) = 0. (A2) 



2m \dr 2 r 2 J 2 
The solution to this is well known and is given by 



rfne 0) 



^ex P (-r 2 /20(y) m < +1/2 (^) (A3) 



with energy E nl = huj(2n + £ + 3/2) where I = ^/KJ mui and n is the number of radial nodes 
in the wavefunction. 

Transformation of (jAl|) into hyperspherical coordinates using Eqs. 12.21 12.31 12.51 and 
12.71 yields a Schrodinger equation that separates into hyperradial and hyperangular pieces. 
The hyperangular solution is a hyperspherical harmonic that diagonalizes A 2 . The resulting 
hyperradial Schrodinger equation is given by 
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where K = A + 3 (N — 1) /2. Comparing this with Eq. IA2I we see that if we make the 
substitutions £ — > A', n — > x, m — ► M, r -> and r/ n£ (r) -> R {ZN - l ^ 2 F xK (R) the 
single particle radial Schrodinger equation becomes the A" particle hyperradial Schrodinger 
equation. With these replacements the solution is evidently 

(R) = A xK exp (-R1/2C) (£) ^ , (A5(a)) 



where C = 1/yN and x is the number of hyperradial nodes in the A" body system. For 
X = this is the same hyperradial wavefunction written in Eq. 12.151 The total energy is 
given by 

E xK = to(2x + K + ^=fu(2 X + \+™ > \ (A5(b)) 

Now that we have a hyperspherical solution we can compare it with the solution to Eq. 
IA1I written in terms of independent particle coordinates (fi, r 2 , rjy). This equation is 
clearly separable in each coordinate r*j, and its solution is a product of A" single particle 
wavefunctions 

TV 

1=1 

Where / ni ^ (r) is given by Eq. IA3[ oji are the spherical polar angular coordinates of the ith 
particle and yp m (u) is a normal 3D spherical harmonic. The energy for this independent 
particle solution is given by 



E = frw 



3N Y 



+ J2 ( 2n i + 



2 

i=i 



(A6) 



Now that we have seen that Eq. IA1I is separable in both hyperspherical and independent 
particle coordinates, we may compare Eqs. |A5(b)| and IA6l to find that 



N 



2x + A = ^(2^ + 4). (A7) 



We should note that this does not mean that the independent particle solution and the hy- 
perspherical solution are the exact same, only that the hyperspherical solution F (R) <3>a (^) 
must be a linear combination of independent particle solutions of the same energy. 



APPENDIX B: PROOF OF THEOREM □ 



Here we will use the results from Appendix A to prove theorem ^ We proceed by 
assuming that there exists a fully antisymmetric, ground state, hyperspherical solution to 
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the Schrodinger equation for N non-interacting fermions in an isotropic trap with energy 
given by Eq. |A5(b)| with % > 0- If we can show that this leads to a contradiction then we 
have, from equation IA71 that there is only one A for all of the ground state configurations. 
From Appendix A we know that this hyperspherical solution must be a linear combination 
of antisymmetric, degenerate, ground state solutions in independent particle coordinates 

F xK (R) $ A (Cl, ax, a 2 , ...a N ) = ^ D u (f h f 2 , ...,f N , a x , a 2 , ...a N ) 

V 

where each D v is of the form of Eq. 12.131 v is used to distinguish between degenerate states 



of the same energy and F x x (R) is given by Eq. A5(a) Again, for brevity, we suppress the 
spin coordinates (01, a 2 , ...ojv) in $a- We now note that the hyperradius R is completely 
symmetric under all transpositions of particle coordinates, thus all of the transposition 
symmetry must be contained in the function $^ (Cl)- From $^ (Cl) we construct a new 
completely antisymmetrized hyperspherical solution Fqk (R) $a (Cl) which has energy 



Use of Eqs. |A5(b)| andlATI gives 



E = !>.c { A + ^ 



where E^j = fiw 



3N/2 + Eti (2n< + 



is the ground state energy as defined by any of 



the functions D v y{ri}f =1 j . Thus our new function F 0K (R) $ A (Cl) has energy 

E — E NI — 2 X 

which would lie below the ground state energy, a contradiction unless x = 0- 

In the above analysis we assumed a degenerate set of solutions at the ground state energy. 
For non-degenerate solutions the proof becomes trivial, as any nondegenerate solution must 
be the same no matter what coordinate system it is expressed in. From this the rest of the 
proof follows in the same way, and the final, unique A for this system is given by Eq. IB H 



A = ^-H. (B2 ) 
hw 2 K ' 
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APPENDIX C: CALCULATION OF THE S-WAVE INTERACTION MATRIX 
ELEMENT 



To calculate CV it is useful to start with a more general interaction. We assume that 
the interaction term in the total iV body Hamiltonian is such that at a fixed hyperradius 
Uint (fij) is separable into a hyperradial function times a hyperangular integral, i.e. 

u int fa) = v ijR (R) v ijn (n) . (ci) 

From properties of the ^-function and Eq. 12.41 it is easy to se that the Ui n t (r) = gS 3 (r) fits 
this criteria. While Vij n (Q) might have some very complex form, it will be seen shortly that 
only the form of V%j R (R) and Ui nt (fv,) will matter. 



U, tJ (R) = ( i, 



5> 

i>j 



int y^ij 



we may again use the antisymmetry of $ A to exchange r*j — > r 2 and r} — > r\ and arrive at 

N(N - 1) 

U eff (R) = 1 2 ' <$ A \U mt (f 21 )| $ a) • 
From Eq. ED we can write that 

U eff (R) = (V R (R) (C2) 

Where ( = ($ A \V n $ A ) is the analogue of gC N /N 3 / 2 in Eq. B~3 To find C we may 
substitute in the definition $ A (^) from Eq. I2.1(i| multiplying by R 3N ~ 1 G (R) 2 on both sides, 
integrating over R and using Eq. IC1I to replace V 2 \ R (R) V 2 \ Q (fi) gives a simple equation 
that may be solve for (. 

(a =f3 (C3(a)) 
a = J R 3N ~ l G {Rf V 2 i R (R) dR (C3(b)) 

N 

/? = iV(i !~ 1) / n ^ (^}f=0 ( c3 ( c )) 

J 3=1 

xU int (^^({r;}^) 

where .D ^{rjj-^J is the Slater determinant defined in Eq. 12.131 and G (R) is defined as in 
Eq. 12.151 We have also used that fact that R 3N ~ 1 dRdQ = Y[ d 3 rj Q|. £ can now be found 

3=1 

as a ratio of two integrals. 
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The integral, a, on the LHS of Eq. |C3(a)1 can be calculated directly. The integral, /?, 
on the RHS of Eq. |C3(a)| is now a diagonal, determinantal matrix element. (3 may be 
drastically simplified by using the orthogonality of the single particle basis functions (for 
details see §6-1). 



TV 

! . 



\J2 [ d 3 nd 3 r 2 (C4) 

x (ri) ^ (r 2 ) C/ int (f 21 ) ^ (fx) ^ (f 2 ) 
- 8m Si m Sj ip* (r 2 ) ip* (ri) C/^ (f 2 i) (n) ^ (r 2 ) 

Where ^ (f) is the ith single particle spatial wave function that appears in the product in 
Eq. I2~T31 i.e. 

Hfc (r) = ^ exp (-r 2 /2Z 2 ) (r/Q* +1 ^ +1/a [(r//) 2 ] y AfB( («) . (C5) 

With this result we can now specify to the s-wave interaction. Following Eq. IC1I we 
may identify U int (f 2 i) < — > gS 3 (r 2 i) and V 2 \ R (R) < — > 1/R 3 . With this and Eq. 12. 15^ the 
integral in |C3(b)| is found to be 



" = / 3^ ■ (C6) 

Evaluating Eq. IC4I begins by integrating the 5-f unction over r 2 . This is simple and we 
can clearly see that the two terms in the integral in Eq. IC4I have a common factor of 
\if)i (ri)| 2 \if}j (r*i)| 2 . Factoring this out gives 

N . 

P = f E i 1 ~ / ^ (n)i 2 Itfc (ri)\ 2 • (C7) 

This sum runs over all spatial and spin quantum numbers in a filled energy shell, we may 
break this up into to factors, a spin sum and a space sum 

1/2 

^=9 E (l-<Wm., 



2 / ^ y -'"sj"^ 

m Si ,m s -=— 1/2 



(ri)| 2 



x / rf3ri (ei^)i 2 J (ei^ 

= » / ^ (E (*)l a ) (E (*)l 2 ) • ( C8 ) 
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where the Greek letters v and \x stand for the set of spatial quantum numbers {n, £, m}. The 
spin sum is trivial and is given by J2m^ m,.=-i/2 ~ ^m ai m a ^\ = 2. We may now calculate 
CV for any given filled energy shell using Eq. IC8I and plugging into the relationship 

C N = £. C9 

g a 

This has been done for the first 100 filled shells the results of which are summarized by Fig. 
^ Cn for a selection of shells is given in table □ 

To extract Cn in the limit as N —>■ oo we may examine the expression for (3 in Eq. IC81 
We may note that the sum p(f) = ^ u \ip u (f)| 2 is the definition of the density of a spin 
polarized degenerate fermi gas of non-interacting particles in an isotropic oscillator. In the 
limit where iV — > oo the trap energy of the non-interacting gas dominates the total energy. 
This means that the Thomas- Fermi approximation becomes exact in this limit. Thus we 
may write that 

^ = ^\-w) i 1 -—) (cl0(a)) 

N m , = I <Prp(r) (C10(b)) 



The chemical potential /x may be found by the condition given in Eq. |C10(b)| where N ms is 
the number of particles with the same spin projection m s . The system we are considering is 
an equal spin mixture so that Nf = vVj_ = N/2. Thus we find that // = tko {?>N) 1 ^ . Plugging 
this into Eq. IC8I gives 

2 , 3 „ 



(3 = gJ [p(r)Yd 3 r 
/2» 

y V 3 315 7r 3 / 3 v ; 



We insert this into Eq. IC9I with Eq. IC6I to find that 



c « - yy ^/ V-v <ci2) 

1 A + 



2 

Using Eq. I3.4al for the large N behavior of A, in the limit N — > oo 
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which is the quoted large N behavior in Eq. 13.61 

It should be said that the same formalism that is presented in this paper can be applied to 
this system in center of mass coordinates. This is done by first choosing an appropriate set 
of Jacobi coordinates, afj = (i + 1) fX)}=i ^j'A — r «+i) f° r * = 1> ■ N — l with the center 
of mass vector defined as x cm = YljLi Fj/y/N. Hyperangular coordinates are now used to 
describe the 3iV — 3 degrees of freedom in the Jacobi coordinates where R 2 = ^2^=1 x 2 /N = 
(^2f=i r | — x "crnj /N- The 3N — 4 hyperangles needed are now defined with respect to the 
lengths of the Jacobi vectors in the same way as in Eqs. 12.31 and 12.41 Under this coordinate 
transformation we find the Hamiltonian to be given by 

H = H cm + Hr,si 

Where Hqm is the hamiltonian for the center of mass coordinate x cm and Hrq is an operator 
entirely defined by hyperspherical coordinates, i.e. 

-h 2 l 

h cm = 7^- V L + 2 muj2x c™ 

With this Hamiltonian we make the ansatz that for Hty = E^fr, ^ = \ (x cm ) G (R) $a (^) 
where \ is a wave function describing the center of mass motion, G (R) is a nodeless hy- 
perradial function and $a (^) is the lowest hyperspherical harmonic in the 3N — 4 angular 
coordinates. We are then looking for the matrix element (<3>a {Hr^I where the integral is 
taken over all hyperangles at fixed R. Theorem 1 still applies with the added idea that \ is 
given by the lowest s-wave state of the center of mass in an oscillator. From here the analysis 
presented in this section for trap centered coordinates still holds with the added change that 
the dimension of the hyperradial integral in Eq. |C3(b)| is three dimensions smaller. This 
leads to a factor in the interaction matrix element given by 



C N -* > ,tfW +HN 



2 / V 2 

For smaller iV this factor changes the interaction considerably, but for larger N it quickly 
goes to 1. Thus, in the large N from limit, besides extracting the center of mass "sloshing" 
modes, there is very little difference from the trap center coordinate systems in this Jacobi 
coordinate formalism. 
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VI. APPENDIX D: CALCULATING p (0) 



N J dm 3 (r N )\$ x 



We start from Eq. 13.141 

p (0) = N J dRR 3N ~ l \F {R)\ 2 J dm 3 (f N ) |$ A (fi)| 2 . 

If we can find J dflS 3 (fjy) \$>\ (^)| 2 then we have the solution. From Eq. 12.41 and properties 
of the 5-function we may say that 

R 3 ' 

We multiply this on both sides by the non-interacting hyperradial function i? 37V_1 (G {R)) 2 
and integrate over R. 

N J J 5 3 (f N )\<l> x (Q)\ 2 (G(R)) 2 R 3N - 1 dRdQ = ^ J (G ' ^ R m ~ x dR (Dl) 

Inserting G (R) from Eq. 12.151 to the right hand side is gives 

f (G(R)) 2 m _, r[A-3(JV-l)/2] 

£ y -wnw R dR - * £3 r[A _ 3iV /2] • (D2(a)) 

From the definition of $a in Eq. 12.161 we see that the left hand side of ID1I is a determi- 
nantal matrix element of a single particle operator integrated over all independent particle 
coordinates. 

N f f 6 3 (r N )\<S> x (n)\ 2 (G(R)) 2 R 3N - 1 dRdn = N f ]J d 3 r \d ({r,}^) |V (r N ) 

(D2(b)) 

Where D (^{ri}f = ^j is the Slater determinant wave function defined in Eq. 12.131 Referring 
to the definition of the density p in Eq. 13.131 we see that this is merely the peak density of 
the non- interacting system p^i (0). Thus 

/ 3 r [A - 3JV/2] 
s HmK ' N S / 2 T [A - 3 (N - 1) /2] ' 
The peak density is then given by 

which is what we were seeking to show. 
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